home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / demos / icse / icseinit.ex < prev    next >
Text File  |  1999-09-16  |  1KB  |  30 lines
  1. // icse.init.bas : initialisation et tests pour icse
  2. //**************************************************
  3. //
  4. //
  5. // creation y0 (etat initial)
  6. if exists('y0')=0, y0=ones(1,ny); end;
  7. // creation de b,fy et fu et passage en vecteur ligne
  8. if exists('b')=0, b=ones(1,ny); end;
  9. if exists('fy')=0, fy1=ones(1,ny*ny);
  10.      else fy1=matrix(fy,1,ny*ny) ; end;
  11. if exists('fu')=0, fu1=ones(1,ny*(nuc+nuv));
  12.    else fu1=matrix(fu,1,ny*(nuc+nuv)); end;
  13. //
  14. format('e');
  15. iu(5)=0;    
  16. [xx,nitu]=size(itu);    
  17. [yy,ndtu]=size(dtu);    
  18. if xx+yy>2, error('itu et/ou dtu n'est pas un vecteur ligne'), end;     
  19. u=[]; if nuc>0, u=uc; end;    
  20. if nuv>0, u=[u,uv]; end;
  21. nu=nuc+nuv*(nti+ntf+1);   
  22. if size(u)<>[1,nu],error('dimensions du controle incompatibles'),end;   
  23. clear xx yy;
  24. //    initialisation du common icsez
  25. [nitv,nrtv,ndtv]=fort('icse0',nu,1,'i',t0,2,'d',tf,3,'d',dti,4,'d',..   
  26. dtf,5,'d',ermx,6,'d',iu,7,'i',nuc,8,'i',nuv,9,'i',ilin,10,'i',nti,..    
  27. 11,'i',ntf,12,'i',ny,13,'i',nea,14,'i',itmx,15,'i',nex,16,'i',nob,..    
  28. 17,'i',ntob,18,'i',ntobi,19,'i',nitu,20,'i',ndtu,21,'i','sort',.. 
  29. [1,1],22,'i',[1,1],23,'i',[1,1],24,'i');  
  30.